APS/123-QED 



O 

(N 

< 



o 
o 

u 

^: 
^ : 

^ ■ 

e 

I 

a ■ 

o : 

o . 



> 
o 

1^ 



m 
o 
o 



X 



Two-orbital model explains the higher transition temperature of the single-layer 
Hg-cuprate superconducter compared to that of the La-cuprate superconductor 
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In order to explore the reason why the single-layered cuprates, La2-a:(Sr/Ba)2:Cu04 {Tc — 40K) 
and HgBa2Cu04+i {Tc — 90K), have such a significant difference in Tc, we study a two-orbital 
model that incorporates the d,,2 orbital on top of the d^2_y2 orbital. It is found, with the fluctuation 
exchange approximation, that the d^2 orbital contribution to the Fermi surface, which is stronger 
in the La system, works against d-wave superconductivity, thereby dominating over the effect of the 
Fermi surface shape. The result resolves the long-standing contradiction between the theoretical 
results on Hubbard-type models and the experimental material dependence of Tc in the cuprates. 

PACS numbers: 74.20.-z, 74.62.Bf, 74.72.-h 



The physics of high-Tc superconductivity, despite its 
long history, harbors rich problems which are still open. 
Specifically, given the seminal discovery of the iron-based 
superconductors [3| and their striking material depen- 
dence of TcLS], it should be important as well as intriguing 
to have a fresh look at the cuprates, which still have the 
highest Tc to date, to understand their material depen- 
dence of the Tc- One of the basic problems is the signifi- 
cant difference in Tc within the single-layered materials, 
i.e., La2-a;(Sr/Ba)2,Cu04 with a maximum Tc of about 
40K versus HgBa2Cu04+5 with a Tc — 90K. Phenomeno- 
logically, it has been recognized that the materials with 
Tc ~ lOOK tend to have "round" Fermi surfaces, while 
the Fermi surface of the La system is closer to a square 
shape which impHes a relatively better nestingjl, 

Conventionally, the materials with a rounded Fermi 
surface have been modeled by a single-band model with 
large second (i2(> 0)) and third (i3(< 0)) neighbor hop- 
ping integrals, while the "low- Tc" La system has been 
considered to have smaller t2,ty,. This, however, has 
brought about a contradiction between theories and ex- 
periments. Namely, while some phenomenological[6] and 
t-J model[7|, [1] studies give a tendency consistent with 
the experiments, a number of many-body approaches for 
the Hubbard-type models with realistic values of on-site 
U show suppression of superconductivity for large ^2 > 
and/or < 0, as wc shall indeed confirm below[9|. 

To resolve this discrepancy, here we consider a two- 
orbital model that explicitly incorporates the (1^2 orbital 
on top of the (1^2 _y2 orbital. The former component has 
in fact a significant contribution to the Fermi surface 
in the La system. We shall show that the key param- 
eter that determines Tc is the energy level difference be- 
tween the dj:2-y2 and d^s orbitals, i.e., the weaker the 
contribution to the Fermi surface, the better for d- 
wave superconductivity, where a weaker contribution of 
the dz2 results in a rounded Fermi surface (which in it- 
self is not desirable for superconductivity), but it is the 



"single-orbital nature" that favors a higher Tc dominat- 
ing over the effect of the Fermi surface shape for the La 
system. 

Let us start with a conventional calculation for the 



single-band Hubbard Hamiltonian, H — X^yo- '-uSo-'-jo' 
^Sj^it^4- Here we take the nearest-neighbor hop- 
ping —ti (~ 0.4eV, see table I) to be the unit of en- 
ergy, U = 6, the temperature T = 0.03, and the band 
filling n ~ 0.85 are fixed, while we vary t2 — — ts with 
t2 > 0. Wc then apply the fluctuation exchange approx- 
imation (FLEX) [lOl. [ll| to solve the linearized Eliash- 
berg equation. Tc is the temperature at which the eigen- 
value A of the Eliashberg equation reaches unity, so A 
at a fixed temperature can be used as a measure for the 
strength of the superconducting instability. We show in 
Figd] A as a function of ( 1 1 2 1 + | i 3 1 ) / 1 ^ 1 1( = 2 1 1 2 1 / 1 i 1 1 here) , 
which just confirm that, within the single-band model, A 
(hence Tc) monotonically decreases with increasing |t2| 
and jiaj. A calculation with the dynamical cluster ap- 
proximation (DCA) shows that a negative works de- 
structively against d-wave superconductivityfl^l, and a 
more realistic DCA calculation that considers the oxy- 
gen orbitals for the La and Hg cuprates also indicates 
a similar tendency [l^. As mentioned above, this seems 
to contradict with the experimental results that the ma- 
terials with larger t2 and have actually higher Tc's CI- 



TABLE I: Hopping integrals within the d^2_y2 orbital for the 
single and two orbital models, and AE = E^2_y2 — E^2. 
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-0.471 


-0.456 


t2[eV] 


0.0284 
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FIG. 1: FLEX result for the eigenvalue, A, of the Eliash- 
berg equation for the single-band Hubbard model plotted as 
a function of (|t2| + lisD/l^il, where we take t2 — — ta > for 
U = 6|ti|, T = 0.03|ti|, and the band filling n = 0.85. Fermi 
surfaces are displayed for two cases (indicated by arrows). 



To resolve this, we now introduce the dx'^_y-2-dz2 two- 
orbital model. For the La system, it has long been known 
that a band with a strong character lies rather close 
to the Fermi energy [l4l - [l6| . More recently, it has been 
discussed in Refs.jl. [l7l| that the shape of the Fermi sur- 
face is determined by the energy level of the "axial state" 
consisting of a mixture of Cu d^^-O Pz and Cu 4s orbitals, 
and that the strength of the contribution causes the 
difference in the Fermi surface shape between the La 
and Hg systems. Namely, the dj,2 contribution is large 
in the La system making the Fermi surface closer to a 
square, while the contribution is small in the Hg sys- 
tem making the Fermi surface more rounded. In Fig|21 
we show the present, first-principles^lS] result for band 
structures in the two-orbital model for the La and Hg 
systems, obtained by constructing maximally localized 
Wannier orbitals [T^. The lattice parameters adopted 
here are experimentally determined ones for the doped 
materials [20l [2l| . We can here confirm that in the La sys- 
tem the main band (usually considered to be the "(i^2_j,2 
band") has in fact a strong d^^ character on the Fermi 
surface near the N point, which corresponds to the wave 
vectors (tt, 0), (0, tt) in the Brillouin zone of the square 
lattice. The d^^ contribution is seen to "push up" the 
van Hove singularity (vHS) of the main band, resulting 
in a seemingly well nested (square shaped) Fermi surface. 
In the Hg system, on the other hand, the d^^ band stays 
well away from Ep, and consequently the vHS is lowered, 
resulting in a rounded Fermi surface. 

If we estimate in the two-orbital model the ratio 
(1^2! + l^aD/l^il within the d^2_y2 orbitals, we get 0.35 
for the La system against 0.41 for Hg (table I), which 
are rather close to each other. This sharply contrasts 
with the situation in which the model is constrained into 
a single band. There, the Wannier orbital has mainly 
dx2-y2 character, but has "tails" with a dz2 character es- 
pecially for the La system. Then the ratio (|i2| + |^3|)/|^i| 
in the single-orbital model reduces to 0.14 for La against 
0.37 for Hg (table I), which is just the conventional view 
mentioned in the introductory part. From this, we can 




FIG. 2: The band structure in the two (d^2_y2-d^2) orbital 
model for La2Cu04 (left) and HgBa2Cu04 (right). The top 
(middle) panels depict the strength of the d^2_y2 (^^2) char- 
acters with thickened lines, while the bottom panels the Fermi 
surfaces (for a total band filling n = 2.85). The inset shows 
the band structure of the three-orbital model (see text) for 
La system, where the the 4s character is indicated. 



confirm that it is the dz2 contribution that makes the 
Fermi surface in the La system square shaped, while the 
"intrinsic" Fermi surface of the high Tc cuprate family is, 
as in the Hg system, rounded. 

Now we come to the superconductivity in the two- 
orbital model. For the electron-electron interactions, it 
is widely accepted that the intraorbital t/ is 7— lOi (with 
t ^ 0.45eV) for the cuprates, so we take U = 3.0eV. The 
Hund's coupling J (= pair-hopping interaction J') is typ- 
ically ~ O.IU, so here we take J = J' = 0.3eV, which 
gives the interorbital U' = U — 2J = 2.4eV. The temper- 
ature is fixed at ksT — O.OleV. As for the band filling 
(number of electrons/site), we concentrate on the total 
n = 2.85, for which the main band has 0.85. Here we ap- 
ply the multiorbital FLEX, as described e.g. in ref . f23| . 
for the three-dimensional lattice taking 32 x 32 x 4 fc- 
point meshes and 1024 Matsubara frequencies. We first 
focus on the La system, and investigate how the dz2 or- 
bital affects superconductivity. Namely, while the on-site 
energy difference, AE = E.j.2_y2 — £^^2, between the two 
orbitals is AE ~ 0.9eV for La2Cu04 (table I), we vary 
the value to probe how the Eliashberg eigenvalue A for 
d-wave superconductivity behaves. The result in Figl3] 
shows that A is small for the original value of AE, but 
rapidly increases with AE, until it saturates for suffi- 
ciently large AE. Hence the superconductivity turns out 
to be enhanced as the dz2 band moves away from the 
main band. Note that this occurs despite the Fermi sur- 
face becoming more rounded with larger AE, namely, 
the effect of the orbital character (smaller dz2 contribu- 
tion) dominates over the Fermi surface shape effect. Con- 
versely, the strong orbital character in the Fermi sur- 
face around the (tt, 0), (0,7r) works destructively against 
(i-wave superconductivity. Physically, the reason for this 
may be explained as follows. First, although the La sys- 
tem has a better nested Fermi surface, we find that the 
strength of the antiferromagnetic spin fluctuations (the 
spin susceptibility obtained in FLEX) in La is only as 
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large as that for Hg. This is intuitively understandable, 
since the two electrons on nearest-neighbor sites are less 
constrained to have antiparallel spins in order to gain ki- 
netic energy when two orbitals are active as in La. Sec- 
ondly, d-wave pairing has a rough tendency for higher 
Tc in bands that are nearly half filled, whereas the 
orbital here is nearly full filled. 

We now focus on how the lattice structure affects AE 
and hence superconductivity. This is motivated by the 
fact that AE should be controlled by the ligand field, 
hence by the height, ho, of the apical oxygen above the 
Cu02 plane [III . To single out this effect, let us examine 
the two-orbital model for which we increase ho from its 
original value 2.4lA with other lattice parameters fixed. 
In Fig|4l[a) , which plots the eigenvalue of the Eliashberg 
equation as a function of ho, we can see that A mono- 
tonically increases with the height. As seen from the 
inset of Fig|4ljb), AE is positively correlated with ho 
as expected, and FigUJb) confirms that the increase in 
A is due to the increase in Ai?[2J] . In these figures, 
we have also plotted the values corresponding to the Hg 
system obtained with the actual lattice structure. We 
can see that, while ho — 2.8A for Hg is larger than 
ho — 2.4A for La, AE ~ 2.2eV for Hg is even larger than 
AE ~ I.3eV, which is the value the La system would 
take for ho = 2.8A. Consequently, A for Hg is somewhat 
larger than that for the La system with the same value of 
ho- This implies that there are some effects other than 
the apical oxygen height that also enhance AE in the Hg 
system, thereby further favoring d-wave superconductiv- 
ity. In this context, the present result reminds us of the 
so-called "Maekawa's plot" , where a positive correlation 
between Tc and the level of the apical oxygen pz hole was 
observed [25}. Since a higher pz hole level (i.e., a lower pz 
electron level) is likely to lower E^^ , the positive correla- 
tion between AE and Tc found here is indeed consistent 
with Maekawa's plot. It can be considered that in La 
cuprates, a considerable portion of the doped holes go 
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into the apical oxygen pz, and this effect is effectively 
taken into account in our model. A more detailed study 
on these issues is now under way, and will be discussed 
in a separate publication. 

Finally, let us discuss the effect of Cu 4s orbital, which 
is the main component of the "axial state" discussed in 
Refs. [1, [l3| . In the present two-orbital model the 4s or- 



bital is effectively incorporated in both of the d^. 



and 



dz2 orbitals, i.e., the Wannier orbitals have tails that have 
the 4s character. In order to make the examination more 
direct, we now consider a three-orbital model that ex- 
plicitly considers the 4s orbital. The band dispersion for 
the La system shown in the inset of Fig|2] shows that 
the 4s band lies well (~ 7 eV) above the Fermi level. 
Nonetheless, the 4s orbital gives an important contribu- 
tion to the Fermi surface in that the ratio (|t2| + l^sD/liil 
within the d^2_y2 sector in the three-orbital model takes 
a much smaller value of 0.10, which should imply that 
it is the path d^^-y-i — >■ 4s — >■ dx2-y2 that gives the ef- 
fectively large t2, is, and hence the round Fermi surface, 
as pointed out previously [3, In this context, it is 

worth mentioning that the path dj.2_y2 ^ dz^ d^^-y^ 
also contributes to t2, ts, but has an opposite sign to the 
4s contribution because the dz^ level lies below dj.2_y2, 
while 4s above fi3.2_j,2[2^. So the two contributions to 
the main band cancel with each other, where the cancel- 
lation should be strong when the energy of the dz2 orbital 
is high as in La. 

We now apply FLEX to the three-orbital model vary- 
ing AE = E,j.2_y2 — Ez2 as in the two-orbital model, 
where we fix the on-site energy difference E^s — Ez2 at 
its original value. We have chosen this because a simi- 
lar three-orbital model constructed for Hg (not shown) 
shows that the on-site energy difference between the 4s 
and dx2_y2 orbitals is smaller than in the La system by 
about leV, so in the Hg system, both of Ex2_y2 — Ez2 
and E/^a — E.j.2_y2 are smaller by about leV, which means 
that the ^^2 and 4s levels shift roughly in parallel relative 
to d^2_y2. It can be seen in Fig 13| that the AE depen- 
dence of A in the three-orbital model resembles that of 
the two-orbital model in the realistic AE range. (When 
AE becomes unrealistically large, i.e., when 4s level is 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 



(a) 



Hg 



(b) 



.0. 

2.4 



2.5 2.6 2.7 

ho[A] 



2.3 2.4 2.5 2.6 2. 7 2. 8 2. 9 0. 8 1 . 
ho [A] 



1.2 1.4 1.6 1.8 2.0 2. 2 

AE(ho)[eV] 



FIG. 3: The eigenvalue. A, of the Eliashberg equation for d- 
wave superconductivity is plotted against AE = E^2_y2 — 
E^2 for the two-orbital (red circles) or three-orbital (purple 
triangles) models for La2Cu04. Corresponding eigenvalues 
for HgBa2Cu04 are also indicated. 



FIG. 4: The eigenvalue of the Eliashberg equation A (red 
circles) when /10(a) or Ai5(b) is varied hypothetically in the 
lattice structure of La2Cu04. Blue diamond indicates the 
eigenvalue of HgBa2Cu04. Inset in (b) shows the relation 
between ho and AE. 
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too close to the Fermi level, the Fermi surface becomes 
too deformed for superconductivity to be retained.) We 
have also calculated the eigenvalue for the Hg system in 
the three-orbital model, and obtained a value very sim- 
ilar to that obtained in the two-orbital model, as plot- 
ted in FiglH If we summarize the three-orbital results, 
while the 4s orbital has an important effect on the shape 
of the Fermi surface, this can be effectively included in 
the dx2-y2 and 1^22 Wannier orbitals in the two-orbital 
model as far as the FLEX studies are concerned. This 
contrasts with the case of the orbital, which, if ef- 
fectively included in the dx2_y2 Wannier orbital to con- 
struct a single-orbital model, would result in a different 
result. This conclusion is natural, since the energy differ- 
ence (~ 1 eV) between d^^-y^ and d^^ orbitals in the La 
system is smaller than the electron-electron interaction, 
which is why the orbital has to be explicitly consid- 
ered in a many-body analysis, while the energy difference 
(~ 7 eV) between d.j.2_y2 and 4s orbitals is much larger 
than the electron-electron interaction, so that the 4s or- 
bital can effectively be integrated out before the many- 
body analysis. So the message here is that the two-orbital 
{dx2-y2-dz2) model suffices to discuss the material depen- 
dence of the Tc in the cuprates. Whether the effect of 
the dz2 orbital can be further incorporated in the on-site 
U or off-site V values (i.e., material-dependent interac- 
tion values) in an effective, single-band model is a future 
problem. 

To summarize, we have introduced a two-orbital model 
to understand the material dependence of Tc in the 



cuprates. We have shown that the key parameter is the 
energy difference between the d^2_y2 and dz2 orbitals, 
where the smaller the contribution of the dz2 orbital, 
the better for d-wave superconductivity, with the orbital- 
character effect superseding the effect of the Fermi sur- 
face shape. It is intriguing to note that the two high 
Tc families, cuprates and iron pnictides, exhibit material 
dependence of Tc that, according to the present study 
and Ref. jl^l , owes to the material dependent multiorbital 
band structures. 

In the present view, the Hg cuprate is "ideal" in that 
the dz2 band lies far below the Fermi level. Neverthe- 
less, there is still room for improvement: as mentioned 
in the outset, within single-orbital systems higher Tc can 
be obtained for smaller ^2 and ^3. It may be difficult 
to make ^2 and ^3 smaller in the cuprates, since they 
are intrinsically large as far as the Cu 4s orbital is ef- 
fective. Conversely, we can predict that materials with 
an isolated single band that has smaller ^2 and ^3 should 
accommodate even higher Tc than the Hg cuprate, pro- 
vided that the electron interaction is similar to those in 
the cuprates. 
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